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POPULATION STUDIES ON EDMONTON MOSQUITOES* 


YOSHITO WADA 

Department of Medical Zoology Quaestiones entomologicae 

Nagasaki University School of Medicine 1:187—222. 1965 

The seasonal fluctuations of each instar larvae and pupae of Culiseta inornata (Williston) 
in a particular pool near the University of Alberta were investigated and an attempt to estimate 
the mortality of the aquatic stages was made. The data for the collections of adults and larvae 
of 26 species of mosquitoes found around Edmonton indicate that the black-legged mosquitoes of 
subgenus Ochlerotatus, genus Aedes are earlier-appearing species than others. The distribution 
pattern of mosquito larvae was firstly demonstrated to follow a negative binomial distribution with 
a common value of constant k for various density levels. Based on this distribution pattern, a 
sequential sampling technique was applied to classify a mosquito population into one of three pre¬ 
defined density levels. This was considered useful in deciding whether or not control is neces¬ 
sary, and in evaluating whether or not control has been successful over a wide area in a relatively 
short time. 


GENERAL INTRODUCTION 

The City of Edmonton has been engaged in the control of 
mosquitoes and has reduced the mosquito population greatly in the 
city (see Klassen and Hocking, 1963 and 1964). However, there are 
still some problems to be solved. They include precisely when and 
how the insecticidal applications should be made for the effective 
and economical control of mosquitoes, how far the larvicide should 
be applied beyond the city limits, and so on. For the settlement of 
them, extensive fundamental studies are required. This report 
deals with the studies conducted in 1964 to approach the problems 
from an ecological point of view. 


BIONOMICS OF EDMONTON MOSQUITOES 

Mosquito Surveys and Identification 

Three types of mosquito surveys were made in 1964. Firstly 
larval (and pupal) surveys were made at pools in various environ¬ 
ments around Edmonton, mostly westward, from April to July. The 
number of dips at each pool was not recorded, except for a few pools 
for determining the distribution pattern of larvae per dip, which will 
be mentioned later. However, care was taken in catching mosquitoes 
so as to represent the mosquito fauna there; only a few dips were 
made at pools with high mosquito density and many dips, sometimes 
more than 50. at pools with low density. 

* Contribution from the Research Institute of Endemics, Nagasaki 
University No. 473 and Contribution No. 142 from the Department 
of Medical Zoology, Nagasaki University School of Medicine. 
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Secondly collections were made of adult mosquitoes, which 
came to feed on me, around a particular pool near the University of 
Alberta at approximately one week intervals. 

Thirdly larval surveys were made at the pool mentioned 
above. The pool harbored almost exclusively Culiseta inornata 
(Williston) and the seasonal changes of immature stages were 
studied. 

The larvae collected were reared in the laboratory to the 
fourth instar or to adults, and identified. Some specimens were 
separately reared to obtain the adults with associated larval skins 
to facilitate determining the species. 

The identification of larvae followed Carpenter and La Casse 
(1955) and Rempel (1950). Adults were identified mostly after 
Carpenter and LaCasse (1955) and Rempel (1953). However, it was 
often difficult to separate them to species, especially rubbed speci¬ 
mens of black-legged female Aedes . In such cases, and even for good 
specimens, the post-coxal scale patch (between the anterior coxa 
and the sternopleuron), mesepimeral scale patch, scales of 
probasisternum, and tarsal claws were useful characters (Beckel, 
1954; Vockeroth, 1954). 

Notes on Some Species 

Aedes communis (DeGeer) and Aedes intrudens Dyar 

A. communis and A. intrudens are black-legged species lacking 
the post-coxal scale patch. The adult female of A. communis is 
usually separable by the contrasting stripes on the scutum from 
A. intrudens with a uniformly colored scutum. However, in some 
specimens of A. intrudens the scutum shows indications of pairedmedian 
brown stripes, and those specimens, particularly when the scales 
on the scutum are not complete, are sometimes hard to distinguish 
from A. communis . 

After examinations of 38 females of A. intrudens and 32 of 
A. communis , some of which were associated with their larval skins, 
it was found that, as described by Carpenter and LaCasse (1955), 
mesepimeral scales reach near lower margin in A. communis , but in 
A. intrudens the lower third or fourth is bare. This seems to be a most 
useful character to separate them. Other characters, which might 
be used, are the number of lower mesepimeral bristles and the color 
of the base of the costa. The lower mesepimeral bristles vary in 
number in both species, but, in the present specimens A. intrudens 
has a smaller number of bristles, ranging from 0 to 3, than A. communis , 
which has 2 to 7 bristles. White scales at base of the costa are 
absent, or if present very few in number, in A. intrudens ; they are 
present in A. communis , 

Aedes hexodontus Dyar and Aedes punctor (Kirby) 

The adults of these two species are very similar to each 
other, however the larvae are distinct. According to Beckel (1954) 
the probasisternum has white scales and an extensive patch of white 
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scales is seen at the base of the costa in A. hexodontus taken in the field 
at Churchill, Manitoba; on the other hand in A. punctor taken there 
scales on the proba si sternum are reduced to a few and there are no 
white scales at the base of the costa or rarely one or two. These 
characters were found useful to separate specimens of these species 
taken near Edmonton also, by examination of females associated 
with their larval skins. 

Knight (1951) recognized two varieties in each species: "type 
hexodontus " and "tundra" variety in A. hexodontus and "type punctor " and 
"tundra" variety in A. punctor . The scutum of females has a broad 
median dark stripe which may be narrowly divided in "type hexodontus " 
and "type punctor", on the other hand in "tundra" variety of both species 
the median dark band is absent or not well defined. 

Of the females of Edmonton hexodontus collected or reared from 
larvae, 9 are "tundra" variety and one is "type hexodontus " variety. 
The latter was collected as an adult on June 7, 1964. In addition to 
these, I have another female specimen of "type hexodontus u variety, 
which was reared from a larva taken near Jasper, Alberta, on May 
16, 1964. The associated larval skin shows that head hairs 5 and 6 
are both double, which agrees with the description given by Knight 
(1951) for "type hexodontus u variety. 

As for A. punctor 9 themany larvae and 18 females, which were 
collected as adults or reared from larvae, are considered all "type 
punctor " variety. 

Aedes niphadopsis Dyar and Knab 

A larva of this species was taken from a collection of small 
scattered pools in a pasture near a creek, about 20 miles west of 
Edmonton, on June 7, 1964, and reared to a female adult. This 

record is new to Canada (Pucat, 1964). 

Aedes pullatus (Coquillett) 

This is a species that lacks the post-coxal scale patch, and 
bears a distinct hypostigial scale patch. The distribution in Alberta 
seems to be limited mostly to mountainous regions. I collected 
many larvae from snow-melting pools in Jasper National Park on 
June 21, 1964, but no specimens were encountered around Edmonton. 

Seasonal Fluctuation and Mortality of Immature Stages of Culiseta inomata (Williston) 

Observations were made on the changes in abundance of each 
instar larvae and pupae of Culiseta inomata (Williston) throughout a 
season at approximately one week interval in 1964 at a pool, ca. 10 
x 3 m, near the University of Alberta. The pool is situated on the 
south bank of the North Saskatchewan river, and receives little sun¬ 
light because of tall vegetation such as poplars around it. For this 
reason, ice remained at the bottom of the pool as late as May 8, 
and the water temperature was relatively low throughout the summer; 
the maximum water temperature was only 18.3 C, on August 17. 

On each day, larvae and pupae were sampled with a dipper 
usually ten times, but when necessary, 20 or 50 times, and the 


190 


Population Studies 


numbers of each instar larvae and pupae were recorded. The popu¬ 
lation of mosquitoes in the pool consisted of only C. inomata , as far 
as the fourth instar larvae were examined. However, from some 
egg rafts collected at the pool on July 6, there emerged some adults 
of Culiseta alaskaensis (^Ludlow) in addition to C. inomata; this indicated 
that a small number of egg rafts, probably one, of the former species 
was mixed in the collection of the egg rafts. Therefore, some 
C. alaskaensis may have bred also in the pool, even so, the number 
seems to have been negligibly small. 

Egg rafts were first encountered on May 25, and oviposition 
continued until August 10. The number of egg rafts per dip and the 
observation for the rafts on the water surface of the pool show that 
the peak of oviposition activity of C. inomata was in the first half of 
June. 

The seasonal distribution for each instar larvae and pupae 
is shown in Fig. 1, The first individuals of larvae in the first, the 
second, the third, and the fourth instar, and pupae were encountered 
on May 25, June 2, June 16, and June 22, respectively. The peak 
in numbers of first instar larvae was June 8, and with the progress 
of the development the time of each peak became successively later; 
the peak for pupae was on July 6. The period between the peaks of 
first instar larvae and of pupae is about one month. This seems to 
be the time required forC. inomata to develop from the first instar 
larva to the pupa; the mean water temperature was 11 C during the 
period. 

The emergence of adults is thoughtto have occurred most 
actively shortly after the peak of pupae, that is in the middle of July. 
This time of peak emergence was ascertained by the fact that many 
pupal skins were observed on the water surface on July 14 and 23. 

It has been reported that the duration of the larval stage of 
mosquitoes such as Anopheles quadrimaculatus Say and Aedes aegypti ( L. ) is 
affected by temperature, nature and amount of food, and density of 
a population (e. g. see Horsfall, 1955). Therefore, the above period 
of one month at mean water temperature of 11 C will be changed to 
some extent according to the conditions in a pool, even when the 
temperature is the same. Also, the remarkable difference in water 
temperatures within a pool (Haufe, 1957) may influence the data. 
However, the difference does not seem to be great, as most larvae 
inhabit similar environments. 

The area surrounded by the abscissa and the curve for each 
instar larvae and pupae in Fig. 1 is dependent on the relative abun¬ 
dance and also on the duration of each instar. In the laboratory at ca. 
23 C, an egg raft of C. inomata was reared to adults, andmean periods 
for each instar larvae and pupae were obtained. If it is supposed that 
these mean periods are kept unchanged also in the present field data, 
we can get the relative abundance by dividing the calculated area from 
Fig. 1 by the mean period. The results are given in Table 1. It is 
recognized from the table that the reduction in the relative abundance 
is remarkable between the first and the second instar larvae, and 
between the third and the fourth instar larvae. The survival rate 
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Fig. 1. Mean number of each instar larvae and pupae of Culiseta momata 
per dip. 
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from the first instar larvae to the pupae is estimated at 63/248 x 
100 = 25%. Thus we get a mortality of 75% for the aquatic stages of 
C. inomata , or slightly higher, as the mortality in the earlier half of 
the pupae is not included in the above calculation. 


TABLE 1 - Relative abundance of each instar larvae and pupae of 
C. inomata in the field. 





Larvae 





1st 

2nd 

3rd 

4th 

Pupae 

Area (days x no. of 

individuals) in Fig. 1 

(A) 

718 

275 

311 

292 

220 

Mean period (days) in 
the laboratory 

(B) 

2.9 

2. 0 

2. 3 

4. 2 

3. 5 

Relative abundance in 
the field 

(A/B) 

248 

138 

135 

70 

63 

The reliability 

of the 

above 

calculation 

depends 

on how 


effectively the material was sampled from the pool and how close 
the relative mean duration for each instar larvae and pupae obtained 
in the laboratory is to that in the field. As will be mentioned later, 
the number of larvae plus pupae of mosquitoes per dip follows a 
negative binomial distribution having a larger variance than a random 
distribution. This means that a larger number of dips is required 
to estimate the population effectively, and the number of dips may 
be too small in the present field data. As mentioned earlier, the 
mean duration of larval stage is affected by temperature, food, and 
population density, but perhaps little affected in pupae by the last 
two factors. Therefore, it is rather difficult to compare the values 
in the laboratory with those in the field. Another difficulty is that 
the temperature in the field changes daily and seasonally. Never¬ 
theless, the above method of estimating the mortality is of value as 
a first approach to this important subject. 

In any case, it seems that the mortality in the aquatic stages 
of C. inomata is fairly high in the field. The factors responsible for 
this are not known. However,, physiological disorder or a sort of 
disease is supposed, as some dead larvae were found and all attempts 
to find predators in the pool failed. 

Seasonal Occurrence of Edmonton Mosquitoes 

Table 2 gives the number of larvae (and pupae) collected 
around Edmontonand the number of collections in which each species 
was found. Mosquitoes were encountered at 30 pools out of more 
than 60 examined. Since the number of dips varies from pool to 
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pool, the number of larvae shown in the table does not represent exactly 
the relative abundance of each species. However, the main features of 
seasonal appearance are clearly seen. 


TABLE 2 - The total number of larvae and pupae collected around 
Edmonton, and the number of collections (within paren¬ 
theses) in which each species was found. 



April 

May 

June 

July 


Species 

early late 

early late 

early late 

early late 

Total 

Anopheles 






earlei 


58(2) 

4(2) 

8(1) 

70(5) 

Culex 






tars alls 




1(1) 

1(1) 

territans 



KD 

5(1) 3(2) 

9(4) 

Culiseta 






alaskaensis 




1(1) 

1(1) 

inomata 



20(1) 436(2) 

7(2) 158(2) 

621(7) 

morsitans 



2(1) 


2(1) 

Aedes 






campestris 

37(1) 




37(1) 

canadensis 

1(1) 




1(1) 

cataphylla 

19(3) 

3(1) 



22(4) 

cinereus 

KD 

3(2) 4(1) 

2(1) 


10(5) 

communis 

10(3) 34(1) 

16(2) 



60(6) 

dorsalis 


1(1) 

1(1) 2(1) 


4(3) 

excrucians 

5(2) 6(1) 

6(1) 13(1) 

2(1) 


32(6) 

fitchii 

31(2) 65(2) 

3(1) 47(1) 

2(1) 


148(7) 

flavescens 


3(1) 

1(1) 


4(2) 

hexodontus 

3(1) 3(2) 

1(1) 



7(4) 

implicatus 

23(3) 24(3) 

20(1) 1(1) 

3(1) 


71(5) 

increpitus 

1(1) 33(2) 

33(1) 



67(4) 

intrudens 

40(2) 


2(1) 


42(3) 

niphadopsis 



1(1) 


1(1) 

pionips 


1(1) 



KD 

punctor 

42(3) 2(1) 

12(1) 

1(1) 


57(6) 

riparius 

29(2) 27(1) 

3(2) 



59(5) 

spencerii 

6(2) 

10(1) 



16(3) 

vexans 


2(1) 

126(4) 


128(5) 

Total 

209(8) 233(5) 

80(4) 160(2) 

37(1) 569(6) 

20(2) 163(2) 






1471(30) 
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The results of the collections of female mosquitoes, which 
came to feed on me, around a pool on the south bank of the North 
Saskatchewan river near the University of Alberta are given in Table 
3. This table also indicates an aspect of seasonal fluctuations of 
mosquitoes. 

From these tables and some other data, seasonal occurrence 
of mosquitoes in 1964 is given below. 


Anopheles 

Anopheles earleiy argas hibernates as an adult female. Many 
larvae were found from late May to early July (Table 2), and one 
female was collected at the campus of the University of Alberta on 
May 26. Most of 58 larvae collected in late May shown in Table 2 
were in the second instar and a few were in the first and a few in the 
third. Thus it seems that hibernated females appear and oviposit 
their eggs from May, and the emergence of adults occurs from June. 
Oviposition continued at least until the beginning of July, as two first 
instar larvae were encountered in early July. 


Culex 

The species of Culex found were C. tarsalis Coquillett and 
C. territans Walker. Both hibernate as adult females. 

Although only one larva of C. tarsalis was collected, the 
hibernated females are considered to oviposit late in the season,as 
it is reported that in irrigated areas of Alberta the larvae are found 
abundantly in July, August, and September (Shemanchuk, 1959), and 
in Saskatchewan the first larvae do not appear until early July 
(Rempel, 1953). 

The larvae of C. territans were collected in late June to late 
July(Table 2), and this seems to be also a late-appearing species . 

Culiseta 

Three species were encountered around Edmonton, namely 
C. alas kaensis (Ludlow), C. inornata (Williston), and C. morsitans (Theobald). 
They all hibernate as adult females. 

The first egg raft of C. alaskaensis was found on July 6, as 
mentioned earlier, and one larva was collected in late July (Table 
2). According to Jenkins (1948), overwinteredfemales were common 
from late April to mid-June and all instars of larvae were found 
fromMay lltoJuly 10 in Alaska. Therefore, the larvae may appear 
earlier than July also around Edmonton. 

Table 2 indicates that the larvae of C. inornata were collected 
from early June, and this agrees with the data mentioned earlier. 
The peak of oviposition was found to be in early June and the 
peak emergence occurred in mid-July. The feeding activity seems 
to be limited mainly to the period from late May to early July, as 
judged from the number of females attracted to man (Table 3), and 
this is justified by the time of the peak of oviposition. Those females 
are considered overwintered ones. However, a small number of 
females oviposited as late as August 10 as mentioned earlier. It is 
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not known whether such oviposition was derived from overwintered fe¬ 
males or from newly emerged ones. 


TABLE 3 - The number of female mosquitoes collected ground a pool on 
the south bank of the North Saskatchewan River near the 
University of Alberta, Edmonton. 



May 


June 



July 


Aug. 

Sept. 















Species 

19 25 

2 

8 

16 

22 

6 

13 

23 

30 

10 25 

11 16 

Total 

Culiseta 













inornata 

1 

1 

1 

1 

1 


1 





6 

Aedes 













cataphylla 

1 











1 

cinereus 





1 

1 




1 3 


6 

communis 

1 


1 

2 

2 







6 

excrucians 





1 







1 

fitchii 





12 

2 

7 

1 



1 

23 

hexodontus 

1 


1 


1 







3 

implicates 

4 3 


1 


25 

1 






34 

increpitus 





3 

1 



2 

1 


7 

intrudens 


1 










1 

punctor 



1 




2 



1 


4 

riparius 





1 


2 

1 


1 


5 

stimulans 











1 

1 

vexans 







3 



3 


6 

Total 

5 6 

2 

5 

3 

47 

5 

15 

2 

2 

6 4 

1 1 

104 


* 

One hour collection was made in the afternoon each day, excepting two 
hour collection on June 22. 


Two laryae of C. morsitans were obtained in early June. Rempel 
(1953) reported the adults in July. This species perhaps spends a similar 
life cycle to C. alaskaensis and C. inornate in Alberta. 


Aedes 

All Aedes species recorded here hibernate as the egg stage. 
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Black-legged species belonging to the subgenus Ochlerotatus are 
generally earlier-appearing species than other mosquitoes. The dates 
of the collections of the larvae and adults in those black-legged species 
(from Tables 2 and 3), together with the records of the larvae and adults 
in Saskatchewan by Rempel (1953) and the dates of emergence near 
Edmonton by Klassen and Hocking (1964) are shown in Table 4. 


TABLE 4 - Summary of the occurrence of black-legged Ochlerotatus . 


Aedes Collections (1) of 

(Ochlerotatus) Larvae Adults 

Records (2) of 
Larvae Adults 

Dates of 
emergence (3) 

cataphylla 

early Apr. May 19 

late Apr. early May 

May 14 


-early May 


-June 15 

communis 

early Apr. May 25 

mid-May 

May 30 


-early May -June 22 


-June 7 

hexodontus 

early Apr. May 25 




-early May -June 22 



impiger 


late May 

May 19 

implicatus 

early Apr. May 19 


May 14 


-early June -July 6 


-June 17 

intrudens 

early Apr. June 2 

June 5 



-early June 

-Aug. 18 


niphadopsis 

early June 



pionips 

early May 

as late as 




mid- July 


punctor 

early Apr. June 8 

May 



-early June -Aug. 25 



spencerii 

early Apr. 

common in abundant 




late Apr. by May 10 



(1) Tables 2 and 3; (2) Rempel, 1953; (3) Klassen and Hocking, 1964. 


It is apparent from the table that the larvae of most species appear 
very early in the season. However, A. pionips Dyar is perhaps a slightly 
later species, as indicated by Haufe (1952) and Rempel (1953), and it 
seems in A. intrudens and A. punctor that the hatching from eggs continues until 
later in the season, or the life span of the adults is longer. As for 
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A. niphadops is it is not clear, as only one larva was collected. 

For banded-legged mosquitoes of subgenus Ochlerotatus, similar data are 
also given in Table 5. 


TABLE 5 - Summary of the occurrence of banded-legged Ochlerotatus 


Aedes 

(Ochlerotatu 

Collections (1) of 
s) larvae adults 

Records (2) of 
larvae adults 

Dates (3) of 
emergence 

campestris 

late Apr. 

early Jun. 

May 19 

(4,5) 




canadensis 

late Apr. 

mid-May 

May 30 



-late Jul. 

-Jun. 7 

dorsalis (4) 

late May 

generally late Jul. 



-late Jun. 

early Jul. -Aug. 


excruciansi 5) 

early Apr. Jun. 22 

early May 

May 27 


-early Jun. 


-Jun. 4 

fitchii • 

early Apr. Jun. 22 

May June 

May 27 


-early Jun. -Sept. 11 

-early Jul. 

-Jun. 17 

flavescens( 4) 

late May 

mid-May generally 

Jun. 7 


-early Jun. 

Jun. - Jul. 

-Jun. 17 

increpitus 

early Apr. Jun. 22 

early May late May 



-late May -Aug. 10 

-June 


riparius 

early Apr. Jun. 22 

mid-late 



-early May -Aug. 10 

May 


stimulans 

Sept. 16 

late May late May 

May 30 



-early Jul. 

-Jun. 4 


(1), (2), and (3) see table 4; (4) a second generation may occur;(5) long-lived 
species, occasional specimens may be encountered late in the season. 


Of the species given in the table, A. excrucians Walker, A. fitchii (Feltand 
Young), A. increpitus Dyar, A. riparius Dyar and Knab, and A. stimulans (Walker) are 
considered woodland species and have only one generation a year. The larvae 
appear as early as most black-leggedmosquitoes, but the emergence is delayed 
because of slower development, as indicated by Haufe (1953 and 1956) and as 
recognized by the fact that the black-legged mosquitoes emerged earlier than 
the banded-legged ones, when the larvae from the same pool were reared in 
the laboratory. The females were collected as late as September llin A. fitchii 
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as August 10 in A. increpitus and A. riparius, and as September 16 in 
A. stimulans . These facts seem to indicate that the life span of adults 
of those species is very long, asRempel (1953) stated that occasional 
specimens of A. excrucians may be encountered in mid-summer. 
A. canadensis (Theobald) is also a wood-loving species. The larvae 
appeared as early as other banded-legged species mentioned above. 
Occasionally hatching occurs in the fall in Illinois (Horsfall, 1955). 

Other tabulated species, A. campestris Dyar and Knab, 
A. dorsalis (Meigen), and A. flavescens (Muller), are grassland-lovers, 
and a second generation may occur, when the environment is favor¬ 
able. They seem to be slightly later-appearing species than the 
woodland species. 

/4e</es vems (Meigen), which belongs to the subgenus 
Aedimorphus , is found in the three main ecological zones in Saskat¬ 
chewan, the prairies, aspen grove region, and coniferous forest 
(Rempel, 1953). This species seems to have multiple generations 
when the conditions are favorable. The larvae were collected in 
early May to late June and the adults on July 13 and August 10. It is 
apparently a late-appearing species. 

Black-legged A. (Aedes) cinereus (Meigen) seems to.be rather 
late in appearance, though the first larva was collected inlateApril. 
The adults were collected from June 22 to August 25. 


DISTRIBUTION PATTERN OF MOSQUITO LARVAE 

Introduction 

Populations of animals may be effectively estimated on the 
basis of their distribution pattern, and much has been published on 
this subject with various kinds of animals, among which however 
mosquitoes are not included. In applying the sequential sampling 
technique, which will be described later, and also in comparing the 
population densities at different pools, it is required toestablish the 
nature of the frequency distribution pattern of mosquito larvae (and 
pupae). 

A dipper is usually used for collecting mosquito larvae, and 
is considered a handy and reliable tool. Here, an attempt has been 
made to analyse the distribution pattern of mosquito larvae in their 
habitats by using the number per dip. 

Collections Used for the Determination of the Frequency Distribution 

Table 6 gives the data of collections of mosquito larvae for 
determining the frequency distribution pattern of the numbers per dip. 
Collections numbers 9 to 24 in the table are the same data as used 
for the seasonal fluctuation of C. inomata described earlier. The table 
indicates that the collections were made at various habitats of various 
sizes during the period covering May 25 to September 30, and the 
mosquito species collected were distributed in the genera Anopheles , 
Culex , Culiseta , and Aedes . The habitats included a grassland pool, a 
woodland pool, a collection of scattered small pools, and the marginal 
part of a creek, and the mosquitoes were found at some times as a 
single species, and at others mixed. 
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TABLE 6 - Collections of mosquitoes in immature stages for 
the frequency distribution pattern. 


Collection 
number Date 

Habitat 

No. of 
dips 

Mosquitoes 

collected 

1 

May 30 

Permanent pool 
in open place 

100 

Anoph. earlei ; 
Aedes spp. 

2 

May 30 

Temporary grass¬ 
land pool 

100 

Anoph. earlei 

3 

June 7 

Collection of small 
pools in pasture 

100 

Culiseta spp. ; 
Aedes spp. 

4 

June 21 

Permanent pool 
in open place 

50 

Anoph. earlei 

5 

June 24 

Marginal part of 
a creek 

60 

Anoph. earlei ; 
Culex territans \ 
Culiseta inornata 

6 

June 24 

Same as No. 3 

30 

Culiseta inornata • 
Aedes spp. 

7 

July 29 

Same as No. 3 

100 

Culex territans ; 
Culiseta spp. 

8 

July 29 

Same as No. 5 

40 

Culex spp. ; 
Culiseta inornata 

9 

to 

24 

May 25 
to 

Sept. 30 

Permanent wood¬ 
land pool(see 
page 189) 

10-50 

Culiseta inornata 


The Relation Between Mean and Variance of the Numbers Per Dip 

2 

In Table 7 the mean,(x), variance (s ), and range of the 
numbers of mosquitoes per dip are given. The means vary from 
0.02 to 39.40, and the variances from 0. 02 to 3975. 34. The minimum 
value of the range for most collections is zero, and the maximum 
value is up to 206. These figures indicate a great variability in 
number of mosquito larvae between the pools and also within each 
pool. 

Several mathematical models have been developed to describe 
the distribution pattern of animal counts. When the distribution is 
considered random, a Poisson distribution is often applied. In the 
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Poisson distribution, the probability for a given positive integer x, 

is given by P(x) = e“ m m x /x t .(1) 

where m is the mean. It is a property of the Poisson distribution 
that the variance is equal to the mean, and the expression 2 (x-x)^/k 
gives a good approximation to x^ with (n-1) degrees of freedom, 
where n is a sample size (Andrewartha, 1961). 


TABLE 7 - Mean, variance, and range of the number s of mosquitoes 
per dip, together with x^ - test for significant departure 
from Poisson distribution. 


Collection 

number 

Mean (x) 

2 

Variance (s ) 

Range 

X 2 

1 

1. 61 

5.47 

0 

_ 

11 

336.60** 

2 

0. 56 

0.89 

0 

- 

4 

157.41** 

3 

0.49 

2.76 

0 

- 

15 

557.37** 

4 

0. 02 

0. 02 

0 

- 

1 

49. 00 

5 

0. 20 

0.82 

0 

- 

5 

161.07** 

6 

14. 60 

509.21 

0 

- 

85 

1011.52** 

7 

1. 00 

2.51 

0 

- 

10 

248.49** 

8 

2.43 

45. 53 

0 

- 

35 

730.86** 

9 

0. 02 

0. 02 

0 

- 

1 

49. 00 

10 

15. 00 

561.11 

0 

- 

76 

336.69** 

11 

35.70 

3975.34 

0 

- 

206 

1002.15** 

12 

33. 30 

1552.54 

0 

- 

159 

885.78** 

13 

32.40 

882.04 

0 

- 

82 

244.98** 

14 

28. 60 

867.82 

4 

- 

80 

273.06** 

15 

39.40 

2590.27 

1 

- 

164 

591.66** 

16 

17.40 

703.38 

0 

- 

88 

363.78** 

17 

19. 30 

368.46 

0 

- 

50 

171.81** 

18 

11.40 

212.93 

1 

- 

44 

168.12** 

19 

6. 30 

58.23 

0 

- 

24 

83.16** 

20 

5. 50 

74.28 

0 

- 

26 

121.59** 

21 

1. 60 

11.60 

0 

- 

11 

65.25** 

22 

0.45 

1.52 

0 

- 

4 

64.22** 

23 

0. 70 

2. 34 

0 

- 

5 

30.06** 

24 

0. 20 

0. 18 

0 

- 

1 

8. 10 


x^ with n-1 degrees of freedom is calculated by (n-l)s^/x. For 
further explanation see text. For collection number see table 6. 

❖ ❖Discrepancy from Poisson distribution is significant at 1% level. 
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It is apparent from Table 7 that the variation of variance is 
much greater than that of mean, and the Poisson distribution does 
not seem to fit the data excepting collection numbers 4, 9, and 24. 
To make sure, the values of 2 (x-x)^/x = (n-1 )$^/x were calculated 
as the fifth column of Table 7; these were highly significant except 
for the above three collections. This shows that a random distri¬ 
bution - Poisson distribution - could not be rejected in the number 
of mosquitoes per dip, when the population density was as low as 
0.20, and discrepancy from Poisson became greater with the in¬ 
crease of the mean, an aggregated type of distribution being indicated. 

As stated by Waters (1959), there will be some field counts 
for which x^ test will show no significant departure from either 
Poisson or an aggregated-type distribution such as negative binomial. 
It seems that non-significant values of x^ in collections 4, 9 and 24 
are attributable to the sparsity of the population and consequent low 
expectation of occurrence of mosquitoes in, individual dips. 

Goodness-of-Fif to the Poisson and the Negative Binomial 

Insect counts in the field are often fitted fairly well by a 
negative binomial distribution (Andr ewartha, 1961; Anscombe, 1949; 
Bliss, 1953), which is one of the aggregated-type distributions. The 
frequency distribution of the negative binomial is given by expanding 
the expression (q-p)~*S where q-p = 1, p = m/k, m is mean, and 
k is a positive exponent. As the variance of a negative binomial 
approaches the mean, or the over-dispersion decreases, k-*o° and 
p-**0. Under these conditions it can be shown that the distribution 
converges to that for the Poisson (Fisher et al. , 1943). 

Goodness-of-fit to the Poisson and the negative binomial was 
tested (Tables 8 to 11) for the data with 100 dips, i. e. collection 
numbers 1, 2, 3, and 7. 

Theoretical frequencies for the Poisson were calculated 
successively by the following formulae. The probability of observing 


zero count, P(0), is 

P(0) = e“ m .(2) 

and the probability of observing (x+1), P(x+1), is 

P(x+1) =mP(x)/(x+l),.(3) 


substituting sample mean, x, for population mean, m. The 
theoretical frequency is obtained by multiplying each probability by 
the sample size, 100. 

The formulae to be used for the theoretical values of the 


negative binomial (Bliss, 1953) are: 

P(0) = (l+m/k)- k .(4) 

and 

P(x+1) = (x+k)mP(x)/ (x+1) (k+m).(5) 


The constant k can be computed by a property of the negative binomial 
that the variance, (J^> is equal to (m+m^/k), where m is mean, 
substituting again sample mean and variance, x and s2, for m and 

cf 2 . 

In all of four examples shown in Tables 8 to 11, highly 
significant departure from the Poisson was demonstrated (p< 0. 001) , 
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which indicates that the distributions cannot be considered random. 
On the other hand, those distributions agree well with the negative 
binomial, except for collection 3, in which some discrepancy from 
the negative binomial is apparent. In this case, 15 larvae per dip 
were recorded once, which is a very high count compared with the 
others. This high count contributes larger variance, whichin turn, 
yields rather small value of k responsible for the discrepancy. 
Generally speaking, the frequency distribution of the numbers of 
larvae per dip seems to agree with the negative binomial. The dis¬ 
agreement with the negative binomial in collection 3 may be attri¬ 
butable to sampling error. 


TABLE 8 - Goodness - of - fit of Collection No. 1 to Poisson and 
negative binomial distributions. 


No. of 

larvae 
per dip 

Observed 

(O) 

F requency 

Hypothetical 

Poisson(P) N. Binom. (N) 

(O-P) 2 

p 

(O-N) 2 

N 

0 

48 


20. 0 

44. 0 

39.20 

0. 36 

1 

17 


32. 2 

20. 8 

7. 18 

0. 69 

2 

11 


25.9 

12. 3 

8. 57 

0. 14 

3 

8 


13.9 

7. 7 

2.50 

0. 01 

4 

4 ' 








8 





5 

4 


7.4 

8. 3 

0. 05 

0. 01 

6 

3 ’ 






7 

1 






8 

1 

8 





9 

2 


0.6 

6. 9 

91.27 

0. 18 

11+ 

1 






Total 

100 


100. 0 

100. 0 

148.77 

1. 39 


*P< 0.001; **0.50< P< 0.75 DF 4* 3** 


Fitting the Negative Binomial Distribution with a Common k 

Comparison between the means of two or more distributions 
are more direct and unequivocal if they have the same relative 
dispersion in terms of k, and two approaches to a common k were 
described by Bliss and Owen (1948). The fir st of them is a regression 
moment estimate applicable to the present data. The following 
calculation is based on Bliss and Owen (1958). 
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TABLE 9 - Goodness - of - fit of Collection No. 2 to Poisson and 
negative binomial distributions. 


No. of 
larvae 
per dip 

Observed 

(O) 

Frequency 

Hypothetical 

Poisson(P) N.Binom.(N) 

(O-P) 2 

p 

(O-N) 2 

N 

0 

66 


57. 1 

64. 6 

1. 39 

0. 03 

1 

20 


32. 0 

22. 5 

4. 50 

0.28 

2 

8 


9. 0 

8. 2 

0. 11 

0. 00 

3 

4+ 

4 " 

2 _ 

6 

1.9 

4. 7 

8. 85 

0. 36 

Total 

100 


100. 0 

100. 0 

14. 85 

0.67 


*P< 0. 001; **0.25< P< 0. 50 


DF 2* 1** 


Two statistics, x' and y* are computed from the mean and 
variance of each component distribution: 

x 1 = x^ - s^/n.(6) 

y' = s 2 - x. (7) 

where n is sample size. Their expectations are given exactly by 

E(x') = m 2 . (8) 

E(y') = m 2 /k.(9) 

Thus (y r -x'/k) has zero expectation. For a single sample, we have 
the ratio 

1/ki = y'/x'.(10) 

as an estimate of 1/k. The variance of (y’-x f /k) is given to order 
1/ r? by 

V = 2m^ (m-k)^[k(k-1) - (2k-l)/n - 3/n^]/(n- l)k^ . . . (11) 

The invariance w = 1/V is of the nature of a weight. If calculated by 
replacing m by x, by x 1 , and k by an empirical trial value of k', 
we can obtain an estimate of 1/k, l/k c , by 

l/k c = 2 (wx'y’) / Z (wx'2).(12) 

as the slope of a linear regression of y' on x', the regression line 
being constrained to pass through the origin (x 1 =0, y' = 0). 

Referring back to the data of Table 7, x'and y'were calculated 
by formulae (6) and (7) for each collection, and the relation between 
them is given in Fig. 2, in log scales so as to show the values with 
great variabilities in one chart. 

Assumed that a proportional relation holds between the two, 
that is given by y' = (l/k)x', then the relation is represented by a 
straight line with an inclination of one in the figure in log scales, 
because log y 1 = log (1/k) + log x'. The data of Fig. 2 satisfies the 
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above assumption very well. This indicates that the relation between 
x' and y’ is represented by a regression line passing through the 
origin, and, in turn the underlying frequency distributions are 
suggested to be the negative binomial with a common k. It is 
interesting that the same trend seems to be shown in the regression 
of y 1 on x 1 between collection numbers 1 to 8 for various species of 
mosquitoes and 9 to 24 for C. inornata (see Table 6), because the 
inclination of the regression line gives the estimate of k, which is 
considered an intrinsic property of the population sampled (Fisher 
et al. , 1943). However, it is likely that the value of k is species 
specific, and further studies are required. 

It is known that in some cases k increases somewhat as m 
increases (Anscombe, 1949; Morris, 1954; Bliss and Owen, 1958). 
So, thevaluesof 1/k^ calculated by equation (10) were plotted against 
mean, x, in Fig. 3, which indicates, however, no appreciable 
relationship between the two. In order to know the exact situation, 
however, the number of dips for each collection seems to have not 
always been sufficient, and further investigations are required. 


TABLE 10 - Goodness - of - fit of Collection No. 3 to Poisson and 
negative binomial distributions. 


No. of 

larvae 
per dip 

Observed 

(O) 

Frequency 

Hypothetical 

Poisson(P) N. Binom.(N) 

(O-P) 2 

p 

(O-N) 2 

N 

0 

77 


61. 3 

83. 3 

4. 02 

0.48 

1 

15 


30. 0 

7. 2 

7. 50 

8.45 

2 

3 

4 

1 " 


7. 3 

3. 3 

1.49 

0. 15 

4 

15+ 

2 

1 

4 

1.4 

6.2 

4. 83 

0. 78 

Total 

100 


100. 0 

100. 0 

17. 84 

9.86 


*P< 0. 001; **0. 001< P< 0. 005 


DF 2* 1** 


Now, a common value of k will be estimated. The statistics 
x' and y 1 for each of the distributions have already been obtained. 
The next step is to get an initial trial estimate of a common k, k 1 . 
As x varies excessively among the collections, a suitable equation 
for k 1 is 

k' = g/2 (y'/x') 


(13) 
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where g is the number of collections. Thus we got k 1 = 0.2822. By 
using this value, l/k c , an estimate of 1/k, was obtained by equation 
(12) and as its reciprocal k c = 0.2947, which does not differ somuch 
from the first trial estimate k 1 = 0.2822. Thus we have estimated 
a common value of k at 0. 2947. If k c should differ appreciably from 
its trial value, k 1 , recalculation is necessary byreplacing the initial 
k 1 by k c . 

The required tests for agreement with a single k c may be 
arranged as an analysis of variance: 


Effect of 

DF 

ss 

MS 

F 

Slope, l/k c 

1 


B o 

B2/S2 

Computed 
intercept 
against 0 

1 

C+b2-B|- 

lo 

Vs 2 

Error 

g-3 

[wy'2]-B2 

S 2 



where = 2 2 (wx'y 1 ) / 2 (wx'^l 

[wx*2] = 2 (wx*2) - 2 2 (wx 1 ) / 2 w, 
[wx’y'^D (wx'y')-2 (wx')2 (wy’)/2 w, 
[wy'2] = 2 (wy’2) - C, 


C = 2 2 (wy')/2 w, 

= [wx'y']2 / [wx'2] 
2 2 (-) = (2 (-)) 2 . 


TABLE 11 - Goodness-of-fit of Collection No. 7 to Poisson and 
negative binomial distributions. 


No. of 

larvae 
per dip 

Observed 

(O) 

Frequency 

Hypothetical 

Poisson(P) N. Binom. (N) 

(O-P) 2 

P 

(O-N) 2 

N 

0 

53 


36. 8 

54.4 

7. 13 

0. 04 

1 

21 


36.8 

21.7 

6.78 

0. 02 

2 

3 

16 

4" 


18.4 

10.8 

0. 31 

2.50 

4 

5 

1 

3 : 

5 

7.7 

9.0 

0.95 

1.78 

6 

10+ 

1 

1 _ 

5 

0. 3 

4.1 

5.39 

0.20 

Total 

100 


100. 0 

100. 0 

20.56 

4.54 


DF 3* 2** 

* P< 0. 001; ** 0. 10< P< 0. 25 
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Fig. 2. Relation between two statistics* x' and y' f defined by equations 
(6) and (7). £ : Collection Nos. 1-8; O *• Collection Nos. 
9 - 24. Collection Nos. 4 and 9 are not shown in the figure , 
because x' =0, y' = 0, and also will be excluded in the later 
calculations, because of indeterminate values of y'/x 1 . 


If a single k c is justified, the F-valuein the first row should 
be clearly significant and that in the second rownot significant. The 
calculated values are shown below: 
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Effect of 

DF 

ss 

MS 

F 

Slope, l/k c 

1 

33.7809 

33.7809 

24.4670** 

Computed 
intercept 
against 0 

1 

4.0297 

4.0297 

2.9186 

Error 

19 

26.2327 

1.3807 




• : 


Fig. 3. Relation between mean (x) and estimate of 1/k (1/k^). 
Collection Nos. 1-8; O : Collection Nos. 9 - 24. 
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The results are highly significant for slope and not significant for 
computed intercept against 0, that is a common value of k is justified. 

Consideration of Reasons for a Negative Binomial Distribution 

I have demonstrated that the number of mosquitoes per dip 
follows a negative binomial distribution with a common k. The 
negative binomial is generated by a distribution that is '’contagious' 1 
in the sense that the presence of one individual in a division increases 
the chance of other individuals falling into that division. However, 
as Andrewartha (1961) stated, agreement with the negative binomial 
does not itself permit any inference about the biology of the mos¬ 
quitoes, though a significant discrepancy from the Poisson series 
disproves the hypothesis of random scatter. In fact, according to 
Bliss (1953) the negative binomial may be regarded as being com¬ 
pounded from a number of Poisson series in which the means vary 
in such a way that they are distributed like x^, and furthermore it 
is possible to imagine a number of other models to explain it. 

The present data are not considered the sum of a number of 
Poisson series with different means, and other reasons should be 
sought. 

One of them which might arise is a dipping error, however, 
its effect seems to be of little importance, or at least, the negative 
binomial distribution is not attributable only to it. 

No habitat of mosquitoes in nature is considered so uniform 
that all parts of it are equally attractive to them. Marginal parts 
of a pool are usually preferable to mosquito larvae, and it is a 
common phenomenon that the spatial distribution of the mosquitoes 
is related to water-plants or overgrown vegetation. Thus the hetero¬ 
geneity of the environment seems to be a great reason for the con¬ 
tagious distribution - the negative binomial. In fact, Hocking (1953) 
observed strong aggregation of the larvae oiAedes communis DeGeer , 
due apparently to the effect of sunlight and temperature gradient in 
the pool. 

Another reason to be considered here may be a gregarious 
habit of mosquitoes. Although this has not beenstudied extensively, 
it seems important in the ecology of mosquitoes. It is commonly 
observed in the laboratory that mosquito larvae show some agg¬ 
regated distribution in a tray, in which the environment does not 
appear to differ appreciably. This habit of aggregation differs in 
intensity with species, and, for example, strong aggregation of 
larvae is frequently seen i r\ Aedes aegypti(L,, ) 9 but it is hardly ever 
seenin Anopheles hyrcanus sinensis Wiedeman. The biological meaning of 
this is not clear at the present time, but is interesting in that it 
may be related to the level of optimum density of larvae. At any 
rate, the intrinsic behaviour of mosquitoes may play some role in 
the contagious distribution. 

In short, the heterogeneity of habitat and possibly a sort of 
gregarious behaviour of mosquito larvae are considered to be res¬ 
ponsible for the negative binomial distribution which is characterized 
by a larger variance than mean. 
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SEQUENTIAL SAMPLING TECHNIQUE 


Introduction 

Sequential sampling can be used for classifying a population 
into one of a number of pre-defined density levels, based on the 
accumulated results of each unit sampled. In classifying animal 
populations, it has been applied to the spruce budworm (Morris, 
1954), whitefish, Coregonus clupeaformis (Mitchell) (Oakland, 1950), the 
lodgepole needle miner (Stark, 1952), and an aphid, Myzus persicae 
(Sulzer) (Sylvester and Cox, 1961). However, it has never been 
applied to mosquitoes. 

The great value of this procedure lies in the fact that it in¬ 
volves a flexible sample size in contrast to conventional sampling 
procedures, and it would frequently be possible to determine whether 
or not a mosquito population requires control, or satisfactory control 
has been obtained, with the expenditure of much less time than would 
have been required if the number of sampling units was inflexibly 
fixed (Knight, 1964). Therefore, it would be reasonable to extend 
this technique to the immature stages of mosquitoes. 

The procedure given by Morris (1954) is mainly followed by 
the present application. 

Density Classes 

As mentioned above, the sequential sampling technique is 
used for classifying a population into pre-defined density levels. It 
is desirable that density classes are determined so as to enable us 
to know from these classes whether or not the mosquito density is 
so high that control operations are necessary, or whether a control 
operation has been successful. 

The density classes may be differently set up according to 
the situation in the city or town concerned. Here, I have classified 
density tentatively into three levels indicated by the critical mean 
number of larvae per dip as follows: 

Densit y Mean number of larvae per di p 

Low 0. 1 or less 

Moderate Between 0. 5 and 2. 5 

High 12. 5 or more 

Density class "high" may be regarded as an indication that the 
mosquito density is so high that control is required, or that a control 
operation has influenced the population but little, and "low" may 
indicate that the density is so low that control is not required, or 
that control was satisfactorily done. "Moderate" is the intermediate 
situation between the two. Although the density is not so high control 
may be desirable if it is early in the mosquito season. 

Of course, the necessity of controlling mosquitoes depends 
not only on the mosquito density in each habitat, but also on the 
relative area of the habitat compared with the whole area, as well 
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as the location of those habitats in relation to city or town to be 
protected from mosquitoes. However, it is still true that population 
density must be determinedat eachhabitat before a decision to control 
is taken. 

Acceptance and Rejection Lines 

To apply the sequential sampling technique to the mosquitoes, 
of which number per dip is consider ed tofollowthe negative binomial 
distribution, it is necessary to find a common value of k fitting all 
the data with different levels of mean, and it has been determined as 
0.2947. 

The next step is to set up alternative hypotheses, Hq and H^, 
from the density classes. To distinguish between low and moderate 
densities at a certain probability level, HQand Hi are that the number 
of larvae per dip is 0. 1 or less and 0. 5 or more, respectively; to 
distinguishmoderate and high they are that the number is 2. 5 or less 
and 12. 5 or more. The values of the constants based on the negative 
binomial distribution at the critical densities under these hypotheses 
are shown in Table 12. 


TABLE 12 - Values of the constants at the critical densities under 
the hypotheses of Hq and H^, based on the negative 
binomial distribution. 


Constant 

Low - 

H 0 

D e n s i t 

Moderate 

Hi 

y 

Moderate 

H 0 

- High 
Hi 

Mean = kp 

0. 1 

0. 5 

2. 5 

12. 5 

p = kp/k 

0.3393 

1.6967 

8.4833 

42.4163 

q = 1 + p 

1.3393 

2.6967 

9.4833 

43.4163 

Variance = kpq 

0.1339 

1.3484 

23.7083 

542.7038 


Each pair of hypotheses is accompanied by two possible 
errors: and /3 are the probabilities of rejecting Hq and H^ at the 

respective critical densities,. Here, both oc and /3 were setat.0. 10. 
A rather large value for error probability seems to be suitable for 
rapid mosquito survey, because it reduces the number of dips to be 
taken at each habitat and enables us to decide whether or not control 
is necessary by a quick evaluation of the population density over a 
wide area in a relatively short time. 

Formulae for the acceptance and rejection lines then are: 


d = sn + hQ.(13) 

and 

d = sn + h^.(14) 






Wada 


211 


where d is the cumulative number of larvae in the first n dips. The 


slope of the lines, s, is 

S = k log (qi/qo) / log (p 1 q () /p 0 q 1 ).(15) 

where qQand are the values of qand Pq and p^ are those of punder 
the hypotheses of Hq and (for actual figures see Table 12), and the 
intercepts of the equations (13) and (14) on the d-axis are 

h 0 = log B / log (p^qo/Poll).< 16 ) 

where B = /3 / (1-ot).(17) 

and 

hi = log A / log (Piq 0 /Po<li).( 18 ) 

where A = (1 - fi) / tx. .(19) 


Thus we get the following formulae as acceptance and rejection lines 
for low versus moderate classes, 
d = 0.2267n - 2.4153 

and 

d = 0.2267n + 2.4153, 
and for moderate versus high 
d = 5.0891n - 24.9138 

and 

d = 5. 0891n + 24. 9138, 

as shown in Fig. 4. This graph may be used in the field to determine 
how many dips should be taken at each habitat in order to define the 
density class within the accepted limits of and . It is helpful to 
visualize each pair of lines as enclosing a band from which the plotted 
points must escape before the density class is satisfactorily defined. 

For example, in collection number 1 mentioned earlier (see 
Tables 6 and 7), the first three dips show no larvae. When zero is 
plotted over each number of dips 1, 2, and 3, it is seen that they are 
within both bands of low-moderate and moderate-high. The fourth 
dip yields two larvae and the fifth none, therefore dips 4 and 5 are 
still within these bands. The sixth dip shows three larvae, so 2 4 3 
= 5 is plotted over dip 6. This is shown to have escaped from the 
bands and to have fallen into the moderate zone, so dipping is dis¬ 
continued. Thus collection number 1 is classified into moderate 
density. If the plotted points had escaped into the area above the 
higher band, the density would be classed as high, and if below the 
lower band, the density would be classed as low. 

The Operating Characteristic Curves 

The operating characteristic curves are useful aids in under¬ 
standing how the plan operates. The curve is calculated from 


L (p) =(Ah - l)/(A h - B h ).(20) 

p =[i - (qo/ < ii) h ]/[(Pi' < io / Po ( ii) h " *].Ui) 


where L(p) is the probability of accepting Hq for any possible level 
of the population mean of kp, A and B are taken from equations (19) 
and (17), and h is a "dummy variable” which may be assigned con¬ 
venient values. 

The operating characteristic curve is shown in Fig. 5 by 
plotting L»(p) against population mean, kp. The left-hand curve is 
for low versus moderate density classes. When the mean, kp, is 
0. 1, the probability of accepting Hq (low density class) is 0. 9; 
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Fig. 4. The acceptance and rejection lines. 
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accordingly the probability of accepting H d (moderate density class) 
is 0, 1. When kp = 0. 5 L(p) = 0. 1 for Hq and consequently 0. 9 for 
H d . At these two levels of kp, the probabilities correspond, of 
course, to those previously set for<x and P> . As kp decreases below 
0.1, L(p) for Hq becomes very low. Whenkpisca. 0. 23, the chances 
of accepting Hq and H d are equal. The curve on the right is used in 
the same way for the moderate versus high density classes. The 
overlapping between the two curves is only at negligible probability 
levels. Thus the probability of considering a low density class 
high, or high density class low, is very small. 

The Average Sample Number Curves 

The average sample number curves can be drawn by plotting 
the values for E(n), the mean number of dips which must be taken, 
against kp, the mean number of larvae per dip, as shown in Fig. 6. 
For different values of kp, E(n) is calculated from 

E(n) =[h d - (h 0 -h ± ) L(p)]/(kp- s).(22) 

where hg, h d , L(p), and s are taken from equations (16), (13), (20), 
and (15), respectively. E(n) does not indicate the number of dips 
which must be taken actually at each pool, but its expectation. 

As would be expected, the peaks of the curves in Fig. 6 occur 
where populations are borderline between low and moderate or be¬ 
tween moderate and high, which indicates that relatively more dips 
are required there. 

Applications of the Sequential Sampling Technique in the Field 

In applying the sequential sampling technique in the field it 
is convenient to use tabulations (Table 13) prepared from theaccep- 
tanceand rejection lines, rather than the lines themselves. Dipping 
is continued until the cumulative number falls into one of the density 
classes. It is apparent from the table that at least 11 dips are neces¬ 
sary for the density to be classed into low, and at least six into 
moderate; if the number of larvae in the first dip is 31 or more, 
the density is classified as high without further dips. 

Table 14 gives the results of applications of the sequential 
sampling technique to the data shownin Tables 6 and 7. Itis demon¬ 
strated that the sequential plan can be used to classify the density 
correctly into one of low, moderate, and high density classes. The 
number of dips required for determining the class in various col¬ 
lections ranged from 1 to 20. When the density is high, the required 
number of dips was rather small, as expected from Fig. 6. This 
is of advantage in field work, because it takes much more time to 
count larvae dipped when the density is higher. 

In sampling, the larvae are required to be dipped all over a 
larval habitat. In a large pool, dividing it into a few portions and 
applying the sequential plan at each will facilitate the work. Suggested 
larval survey form is given in Table 15. 

This technique can be used effectively for the evaluation of 
the application of larvicides in a relatively short time. If the control 
operation is successful, then the densities at all pools will fall into 
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Fig. 5. The operating characteristic curves for low versus moderate 
density classes (left) and for moderate versus high (right), 
kp = mean no. of larvae per dip; L (p) = probability of accept¬ 
ing Hq hypothesis. 



Fig. 6. The average sample number curves for low versus moderate 
density classes (left) and for moderate and high (right), 
kp = mean no. of larvae per dip: E (n) = mean no. of dips to be 
taken. 
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the low density level. Also, this may be used for determining whether 
or not a second larvicide application is required specifically for the 
later-appearing mosquitoes. Neces sity for mosquito control depends 
on the productivity of mosquitoes in a particular area, rather than 
the population density at each pool. To approach this, the following 
procedures may be appropriate. Firstly we determine the density 
class at each pool by the sequential sampling technique. Then, we 
take 0, 1, and 10 as indices for low, moderate, and high density 
levels, respectively, and multiply the index by the area of the pool 
(the area of themarginal parts if the larval distribution is confined 
there). If these are summed for a district to be examined, then it 
will represent the productivity of mosquitoes there. The sequential 
plan may be used for comparing regional differences of mosquito 
abundance, which provide us with the knowledge as to which region 
should be stressed for larval control operations. 


TABLE 13 - Sequential table for use by field parties, prepared 
from the acceptance and rejection lines (Fig. 4). 


No. of 
dips 

Cumulative number of larvae 
Low Moderate 

High 

1 





31 

or 

more 

2 





36 

or 

more 

3 





41 

or 

more 

4 





46 

or 

more 

5 





51 

or 

more 

6 


4 

to 

5 

56 

or 

more 

7 


5 

to 

10 

61 

or 

more 

8 


5 

to 

15 

66 

or 

more 

9 


5 

to 

20 

71 

or 

more 

10 


5 

to 

25 

76 

or 

more 

11 

0 

5 

to 

31 

81 

or 

more 

12 

0 

6 

to 

36 

86 

or 

more 

13 

0 

6 

to 

41 

92 

or 

more 

14 

0 

6 

to 

46 

97 

or 

more 

15 

0 

6 

to 

51 

102 

or 

more 

16 

0 to 1 

7 

to 

56 

107 

or 

more 

17 

0 to 1 

7 

to 

61 

112 

or 

more 

18 

0 to 1 

7 

to 

66 

117 

or 

more 

19 

0 to 1 

7 

to 

71 

122 

or 

more 

20 

0 to 2 

7 

to 

76 

127 

or 

more 


Continue to dip until the cumulative number falls into one of the 3 
density classes of low, moderate and high. 
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TABLE 14 - Application of the sequential sampling technique to the 
data shown in Tables 6 and 7 of Section 3. 


Collection 

number 

Mean 
no. of 

larvae 

Density 

class 

determined 

No. of 
dips 

required 

1 

1. 61 

Moderate 

6 

2 

0.56 

Moderate 

20 

3 

0.49 

Low 

11 

4 

0. 02 

Low 

16 

5 

0. 30 

Moderate 

18 

6 

14.60 

High 

9 

7 

1. 00 

Moderate 

20 

8 

2.43 

Moderate 

7 

9 

0. 02 

Low 

11 

10 

15. 00 

High 

1 

11 

35.70 

High 

2 

12 

33.30 

High 

2 

13 

32.40 

High 

2 

14 

28.60 

High 

1 

15 

39.40 

High 

1 

16 

17.40 

High 

7 

17 

19.30 

High 

3 

18 

11.40 

High 

3 

19 

6.30 

High 

2 

20 

5. 50 

Undetermined* 

£ 11 

21 

1.60 

Moderate 

9 

22 

0.45 

Moderate 

6 

23 

0.70 

Moderate 

7 

24 

0.20 

Undetermined* 

£ 11 


During 10 dips made, the density class was not determined. 


SUGGESTED STUDIES TOWARD 
BETTER CONTROL OF EDMONTON MOSQUITOES 


Introduction 

In this section, only ecological questions are discussed, 
although studies are also needed on the identification of mosquitoes 
including the larvae in younger instars, the development of insec¬ 
ticidal resistance, the methods and evaluation of applications of 
chemicals, the effective and economical dosages of larvicides and 
adulticides, the residual effects of insecticides when applied to the 
habitat in the field, and so on. 
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TABLE 15 - Suggested mosquito larval survey form for the appli¬ 
cation of sequential sampling technique in the field. 


MOSQUITO LARVAL SURVEY FORM 

Collection No. Collector: 

Place: Hour: , a.m. p. m. 

Date: , 19 

Breeding place 

-type: woodland - pool, grassland - pool, roadside - ditch, 
small pools in pasture, creek, other ( ) 

- permanent, temporary 

- size 

- notes (marginal vegetation; water plants; animals; 

temperature, pH, cleanness of water; etc.) 


No. of No. of Cumula- 

dips larvae tive no. 

No. of No. of Cumula- 

dips larvae tive no. 

1 

11 

2 

12 

3 

13 

4 

14 

5 

15 

6 

16 

7 

17 

8 

18 

9 

19 

10 

20 


Density class determined: Low, Moderate, High 
Instar of larvae: 

Species identified: 


The Time of Hatching and Emergence 

The prediction of the emergence time of mosquitoes is re¬ 
quired to determine the appropriate time for chemical control. The 
best time for controlling mosquito larvae is before they begin to 
pupate, the pupae being much more resistant to insecticides than 
larvae, but not before hatching is complete. Strictly speaking, the 
above situation is hard to realize in the field, because the time of 
hatching differs between species and also within species so that there 
remain some eggs of Ajedes to be hatched later in the season after 
some adults have emerged. Thus the most effective time for insec¬ 
ticidal applications against mosquito larvae is our special concern. 
For this purpose, many points remain to be studied. These include 
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the studies on the time of oviposition and the durations of egg and 
larval stages. 

In mosquitoes belonging to the genera Anopheles, Culex ,and 
Culiseta, which overwinter as adults, the time of oviposition depends 
on the time of blood feeding and the duration of egg development. 
Blood feeding is certainly related to temperature and possibly to 
adult diapause. The temperature apparently influences the matur¬ 
ation of eggs. 

All Aedes mosquitoes found around Edmonton overwinter as 
eggs. According to Clements (1963), the different Aedes species fall 
fairly clearly into those whose eggs enter diapause and require re¬ 
activation, and those whose eggs merely become quiescent and hatch 
shortly after exposure to an adequate hatching stimulus, although 
they may require a few hours conditioning. Obligatory diapause in 
the egg stage is found in Aedes hexodontus (Beckel, 1958), in Aedes squamiger 
(Telford, 1958), and in Aedes stimulans (Horsfall and Fowler, 1961), 
where exposure to low temperature is required before egg diapause 
can be broken. These mosquitoes have only one generation a year. 
Multivoltine species have facultative diapause, as in Aedes dorsalis 
(Khelevin, 1958), Aedes nigromaculis (Telford, 1963), and Aedes triseriatus 
(Baker, 1935), or have no diapause. 

Most mosquitoes found around Edmonton have one generation 
a year. However, there is a possibility that a second generation 
occurs in some species, such as A. campestris, A. dorsalis , or 
A. flavescens , perhaps in August when the conditions are favorable. 

It is very likely that there is a wide variability in hatching 
response of eggs, sothat the time of hatching has a wide range, even 
for eggs from the same batch. 

Beckel (1958), Telford (1963) and others have discussed the 
mechanism and ecological significance of egg diapause in mosquitoes, 
and much has been published on the hatching response in quiescent 
eggs of Aedes aegypti and some other Aedes species (see Telford, 1963). 
However, the situation is still not clear for most Aedes mosquitoes. 

After hatchingfrom eggs, the development of larvae depends 
on various factors. The most important are temperature, quality 
and quantity of food, and larval density. It is expected that the 
relation between larval period and temperature is described by an 
equilateral hyperbola, or the relation between developmental speed 
and effective temperature (temperature minus developmental zero 
point) is linear, at least within a reasonable temperature range, 
provided that other factors than temperature are constant. Based 
on this relation, Haufe (1953) and Haufe and Burgess (1956) attempted 
to predict dates of emergence in mosquitoes at Fort Churchill, 
Manitoba, and stated: n The tundra species of mosquito ( A. impiger 
and A. nigripes) had lower thresholds of development approximating 
34F; the forest species (A. communis, A.punctor, A. excrucians) hada range 
of 38 - 40 F, excepts, hexodontus . The products of time and temper¬ 
ature for the period of development of both tundra and forest species 
were lower for the smaller than for the larger species". Studies of 
this sort are desirable for all the mosquito species found abundantly 
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around Edmonton. 

It is to be noted that the threshold of development obtained 
from the above relation is slightly higher than the actual value in the 
development of most insects, and is not necessarily the same as the 
critical temperature for hatching. Also, the developmental speed 
differs greatly according to factors such as quality and quantity of 
food and larval density. Thus, for the prediction of the date of 
emergence, careful investigations are required in the laboratory for 
each species. Another aspect to be involved is the relation between 
temperature in pools in various situations and meteorological 
records, for example see Haufe and Burgess (1956) and Haufe(1957). 

Flight Range 

Southwood (1962) stated: "It is suggested that animal move¬ 
ments fall basically into two types: trivial and migratory. Trivial 
movements are normally confined to the territory or habitat of the 
population to which the animal belongs, migratory movements carry 
the animal away from this area. Although there is undoubtedly no 
sharp line but a gradation between these two types, they can be 
distinguished by various ecological, physiological and behavioural 
characteristics", and "The ideal evidence of migratory movement 
is that while engaged in it the animal does not respond to food, a 
mate or habitat, and moves from the actual territory where it has 
developed into an inhospitable terrain: such movement is normally 
at the start of adult life". Provost (1957) reported the findings of a 
mark - and - release experiment with Aedes taeniorhynchus as follows: 
"Migration occurs the night of departure only, therefore twilight 
departures will result in longer migrations than middle-of-the-night 
departures. Appetential (trivialof Southwood, 1962) flights expand 
the range of occupation by a brood much beyond what is established 
by the migration. " Thus mosquito dispersal consists of two phases 
of movement, and its range depends greatly on the migratory flight 
and to a lesser extent on the appetential flight. 

In the appetential flight of mosquitoes, the distribution of 
breeding, resting, feeding, and oviposition sites, and in some 
species overwintering sites, will influence the degree of dispersal, 
because inmosquitoes these sites are situated quite often at different 
places and it is suggested that, within limits, the closer these are 
situated, the shorter the flight range. This should be considered in 
the field data, particularly when mark-and-release experiments are 
conducted. 

As mentioned earlier, 14 species of adult female mosquitoes 
were collected around a pool near the University of Alberta. Ex¬ 
cepting Culiseta inornata , all these mosquitoes are considered to have 
entered from outside or marginal parts of the City of Edmonton, 
since there are no breeding places for these species in the central 
part. This means that they dispersed at least a few miles, and 
should be consider ed potential pests for the Edmonton area according 
to their abundance. Here, it is required to determine the range of 
dispersal for each species. It is not known whether migratory flight 


220 


Population Studies 


was involved in the dispersal or not. It seems reasonable to suppose 
that the range of dispersal is longer inmosquitoes where migratory- 
flight is involved than in others. The investigation of this subject 
will help decide how widely insecticidal application or other control 
of larvae must be extended to control mosquitoes near Edmonton. 

Fluctuation in Numbers 

The habitats of most mosquito larvae are characterized by 
their unstableness. In years with small precipitation, the habitats 
will be greatly reduced in extent, though the amount of standing 
water is influenced by the dryness of the land to some degree, as 
indicated by Rempel (1953), and the reverse is also the case. The 
change of the habitats determines the area of breeding and oviposition 
places available for mosquitoes. Also, larval mortality has aclose 
association with the amount of precipitation in some circumstances, 
since it is often observed that pools dry up befure mosquitoes emerge. 
Thus mosquito abundance is expected to be closely correlated to the 
amount of precipitation. Temperature is also an important factor 
influencing mosquito populations, as indicated by Rempel (1953). 

From the above statement, it is clear that there exists a 
close relationship between mosquito abundance and meteorological 
factors. However, the situation will differ from species to species, 
as evidenced by the fact that some species appear abundantly in one 
year and others in another year. The analysis of these correlations 
over a long period will help in studies on the population dynamics. 

Another approach to studies on population dynamics is the 
analysis of mortality factors in the field, and also the influence of 
various environmental conditions such as trophic factor s, population 
density, and so on, on the fecundity of adults in the field and also in 
the laboratory. 
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